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ABSTRACT 



Aims. Accurate distances to evolved stars with high mass loss rates are needed for studies of many of their fundamental properties. 
However, as these stars are heavily obscured and variable, optical and infrared astrometry is unable to provide enough accuracy. 
Methods. Astrometry using masers in the circumstellar envelopes can be used to overcome this problem. We have observed the OH 
masers of a number of Asymptotic Giant Branch (AGB) stars for approximately 1 year with the Very Long Baseline Array (VLBA). 
We have used the technique of phase referencing with in-beam calibrators to test the improvements this technique can provide to Very 
Long Baseline Interferometry (VLBI) OH maser astrometric observations. 

Results. We have significantly improved the parallax and proper motion measurements of the Mira variable stars U Her, S CrB and 
RR Aql. 

Conclusions. It is shown that both in-beam phase-referencing and a decrease in solar activity during the observations significantly 
improves the accuracy of the astrometric observations. The improved distances to S CrB (418^^g pc) and RR Aql (633^ pc) are fully 
consistent with published P-L relations, but the distance to U Her (266^2^ pc) is significantly smaller We conclude that for sources 
that are bright and have a nearby in-beam calibrator, VLBI OH maser astrometry can be used to determine distances to OH masing 
stars of up to ~ 2 kpc. 

Key words, masers - stars: circumstellar matter - stars: individual (U Her, S CrB, RR Aql) - stars: AGB and post- AGB - techniques: 
interferometric - astrometry 



1. Introduction 

VLBI astrometric observations of circumstellar OH masers 
can be used to determine the proper motion and parallax of 
enshrouded maser bear ing stars (van La ngevelde et al.L llOOdt 
IVlemmings et all 120031, , hereafter vLOO and V03). Obtaining 
accurate distances to these stars allows the inclusion of the 
more extreme Mira stars in studies of the fundamental prop- 
erties of these stars, like the pulsation and mass-loss mecha- 
nism. Accurate parallaxes also allow the growing number of op- 
tical and infrared interferometric observation s of AGB stars to 
' be us e d to determine ac t ual ste llar sizes (e.g. iMennesson et al.L 
120021: IWittkowsM et al.L |2005|) . Currently, most investigations 
of pu lsating stars are dependent on Hipparcos ( Perrvman et al], 
Il997h . which is biased against the highly obscured stars with 
high mass loss. Unfortunately, as these stars are variable, the 
Hipparcos distances are often very uncertain. Maser astrometry 
can remedy this situation. In order to use the maser positions 
to monitor the stellar trajectory, an assumption has to be made 
about the motion of the masers with respect to the star. In V03, it 
was shown that for several stars the brightest, most blue-shifted 
circumstellar maser spot corresponds to t he Amplified Stellar 
Ima ee, as was previously hy pothesized by iNorris et"al] d 19841) 
and lSivagnanam et al.l ( Il990i) . However, not all stars show such 
a maser spot, depending on for instance an inho mogeneous dis- 
tribution of the masing gas dVlemmingsl |2002|) . Still the obser- 
vations in V03 indicate that even without an amplified stellar 
image, VLBI astrometry of OH masers yields highly accurate 
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parallax and proper motion results. The only assumption needed 
for this is that these spots have linear motions in the shell. 

The technique of VLBI astrometry has undergone extensive 
developments in the past several years. The feasibility of annual 
parallax measurements has now been proven at several differ- 
ent frequencies. At frequencies around 1.6 GHz, current astro- 
metric observations provide pulsair velo c ities and distan ces out 
to several kpc (e.g. lBrisken et all [20021: iChatteriee et aUl2005l) . 
At higher frequencies, H2O and methanol maser astrometry is 
being used to obtain star-forming region dist ances up to sev- 
eral kpc with close to ~ 2% accuracy (e.g. Hachisuka et all 
i2006; Xu et all l2006h . H2O maser observations have also re- 
cently been used to mea sure the parallax of the Mira star UX Cyg 
dKuravama et aUl2005h and additional H2O and SiO maser AGB 
astrometric observat ions are expected from the VERA project 
jHonmaet al.Ll2003h . 

In this paper we present a follow-up on our earlier success- 
ful VLBI OH maser astrometry campaign (vLOO, V03). In § |2] 
we discuss the observations, data reduction and error analysis. 
In § [3] we give the results of our parallax and proper motion 
determination. The results are compared with previous distance 
measurements in § |4] along with a discussion about the internal 
maser motions in U Her. The conclusions are presented in § |5] 



2. Observations & Data reduction 

The positions of the 1.6 GHz circumstellar OH masers of a sam- 
ple of 6 AGB stars were monitored at 6 epochs spaced by 2- 
3 months between February 2004 and April 2005. The observa- 
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Fig. 1. The position of the brightest 1665 MHz maser spot of RR Aql with respect to the in-beam calibrator (NVSS J195655- 
013615). The error bars on the positions indicate the position fitting errors described in § 12.41 (left) The solid line is the best fitting 
parallax and proper motion trajectory for the 5 epochs with usable observations, (right) The best fitted parallax signature after 
subtracting the fitted proper motions. 



Table 1. The sample 



Source 


Period 




Hipparcos Parallax 


Hipparcos Proper motion 






(days) 


(kms-') 


(mas) 


RA, dec (mas/yr) 


(pc) 


SCrB 


360 


0.0 


2.40 ± 1.17° 


-8.33 ±0.93 -11.55±0.62 


470 


UHer 


406 


-14.5 


1.88 ± 1.31° 


-16.84 ±0.82 -9.83±0.92 


380 


RR Aql 


394 


27.8 


2.48 ± 2.57 


-24.01 ±4.18 -47.66±2.80 


540 



° recalculated by Knapp et al. ("2003) 
' from Whitelock & Feast (2000) 



tions were perfomed using the NRACQ VLBA and the sample of 
stars consisted of the Mira variable stars U Her, S CrB, RR Aql, 
R Aql, R Cas and Y Cas. Only for the first three of these sources 
did we manage to detect maser emission that was strong enough 
to image accurately. Their periods, velocities with respect to the 
Local Standard of Rest (LSR), Hipparcos parallaxes and proper 
motions and P-L distances are given in Table [T] The observa- 
tions of all sources were scheduled in the same observing block 
for a total of 12 hours. This resulted in approximately 2 hours 
of observations per maser source, including the time spent on 
the extragalactic phase reference sources. Due to some obser- 
vational problems (strong interference, telescope and recording 
failures), a few hours of observations at 3 of the epochs were un- 
usable. Fortunately, none of the three detected sources (U Her, 
S CrB and RR Aql) were affected more than once and we ob- 
tained 5 epochs of good data for each of them. The resulting 
average beam size for the observations was 14 x 8 mas. 



' The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



2.1. In-beam calibration 

To improve upon the results obtained in V03, we used in-beam 
calibration in addition to the r e gular nodding phase referencing. 
As shown in IChatteriee et al.l (l2005h . the astrometric accuracy 
that can be obtained by phase referencing is strongly depen- 
dent on the separation between the target and reference source. 
After initial calibration on a primary (nodding) phase reference 
source, residual calibration errors arise primarily from the un- 
modeled ionosphere between target and reference source. These 
errors can be significantly reduced by using an secondary cali- 
brator in the primary telescope beam, reducing t he angular throw 
as we ll as the need for time extrapolation (e.g. iFomalont et aU 
1999). Such a secondary in-beam calibrator can be much fainter 
than the primary calibrator, with a flux > 10 mJy. In an effort to 
find in-beam calibrators with sufficient flux that additionally are 
suitably compact at VLBA baselines, we observed the bright- 
est NVSS sources dCondon et all 1 19981) near a number of maser 
stars with the Very Large Array (VLA) at 8.4 GHz. In addition 
to the stars making up the sample presented here, we also ob- 
served possible calibrator sources near W Hya and RT Oph. In 
total 40 sources were observed for 3 min each. These were then 
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Fig. 2. Similar to Fig.[T] The position of the brightest 1667 MHz maser spots of S CrB with respect to J 1522+3 144, which was used 
as in-beam calibrator during the 5 epochs presented here. Drawn is the best fitting parallax signature after subtracting the proper 
motions for the 5 new epochs with in-beam calibration (left) and the proper motion and parallax trajectory including the previous 
observations of the brightest 1665 and 1665 MHz OH maser spots (V03) for a total of 13 epochs (right). 



reduced using AIPS without any special processing. The images 
of all the sources, including those that were not detected, are 
presented in Figs.|6ll2l[l]and|9]in the Online material. There we 
also present Table. [3] with positions and fluxes extracted using 
the AIPS JMFIT task. From the VLA observations we deter- 
mined the best VLBA in-beam candidates. No good candidates 
were found for W Hya and RT Oph and also for U Her there was 
no obvious good in-beam candidate. However, as U Her is the 
longest astrometrically monitored OH maser source in was still 
included in our sample even though in-beam calibration turned 
out to be impossible. As in V03, we used J1630-I-2131 with a 
separation of 2.8° as primary calibrator for U Her. For S CrB, 
the strongest in-beam calibrator source was also the primary cal- 
ibrator (J1522-H3144) used in V03 at a distance of 24', so we 
were able to accurately calibrate the masers of S CrB without 
nodding the telescopes. For RR Aql we used both a primary cal- 
ibrator (J2015-0137 at 4.4° from RR Aql) and an in-beam cal- 
ibrator (NVSS J195655-013615), with a target-calibrator sepa- 
ration of ~ 20' (see Table [3] in the on-line material for detailed 
coordinates). 



2.2. Data correlation 

We used a mixed bandwidth setting to observe both the 
1667 MHz OH maser line and the continuum calibrator simul- 
taneously. One frequency band with a bandwidth of 500 kHz 
was centered on the stellar velocity and correlated with moderate 
spectral resolution (1.95 kHz = 0.36 km s '). Simultaneously 
three 4 MHz wide bands were recorded to detect the contin- 
uum reference sources, with one of the bands also covering the 
1665 MHz OH line. These bands were correlated twice with a 
spectral resolution of 15.6 kHz, equaling 2.88 km s The first 



of the wide band correlator passes used the position of the in- 
beam reference source for calibration purposes, while the second 
pass used the stellar position to allow us to detect the 1665 MHz 
OH maser line. 



2.3. Data calibration 

The astro me trie VLBA data was calibrated using AIPS, with 
amplitude calibration based on the antenna system temperatures 
and phase, rate and delay calibration based on the primary and 
in-beam calibrators. We applied ionospheric corrections using 
the AIPS task TECOR, but found that these do not significantly 
improve our data quality. Additionally, we used CLCOR to cor- 
rect for errors in the Earth Orientation Parameters (EOP) that af- 
fected VLBA observations between May 2003 and August 2005 
(Walker et al. 2005). For a consistent solution we determined, 
when possible a model of the primary and available in-beam 
calibrators using all epochs of data. As the primary calibrator 
(J1630+2131) of U Her showed structural changes between the 
different epochs, this forced us to construct independent cali- 
brator models for every epoch. The models were then used to 
iteratively calibrate the data. 

As the calibration solutions from both the primary and the in- 
beam calibrators are determined on the wide-band (3x4 MHz) 
data, special care needed to be taken connecting these solu- 
tions to the narrow-band spectral line data. A special task writ- 
ten in AIPS for the previous observations presented in vLOO 
and V03 was rewritte n in the ParselTongue scripting language 
dKettenis etaTll2006h and was then used to transfer the calibra- 
tion solutions. 
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Fig. 3. Similar to Fig. [T] The position of the most blue-shifted 1665 and 1667 MHz maser spots of U Her, that have been found to 
correspond to the stellar image, with respect to the primary calibrator Jl 630+2131. Drawn is the best fitting parallax signature after 
subtracting the proper motions for combining the results on both maser transitions for the 5 recent epochs of observations (left), and 
the proper motion and parallax trajectory combining the new observations with the 12 previous epochs (right). 



2.4. Error analysis 

For the two sources with in-beam calibrators (S CrB and 
RR Aql), the astrometric errors on the observational data points 
were taken to be the uncertainty in calibrator and maser feature 
position added in quadrature. The uncertainty in calibrator posi- 
tion is determined using Beam/ [2 x SNR^,] for unresolved cal- 
ibrator sources , while that of the maser feature and of the re- 
solved calibrator sources is 6cj„/[2 x SNRc ,,,] with SNRc,,„ the 
signal-to-noise ratio on the calibrator and maser and 0c ,„ the full 
width half-maximum size of the calibrator and the maser feature 
respectively. This size is equal to the Beam-size when the maser 
is unresolved, which is typically the case. The total positional 
error is of the order of 0.2 mas during the best observing condi- 
tions, with the error as a result of the uncertainty in calibrator po- 
sition typically smallest by an order of magnitude. To total error 
increases to ~ 0.7 mas for the weakest maser-calibrator combi- 
nations and the worst ionospheric conditions. The here described 
errorestimates produced least ;^'^-fits with reduced = 1.01 and 
1 .86 for 5 degrees of freedom for RR Aql and S CrB respectively. 

For U Her, where no suitable in-beam calibrator was found, 
astrometric errors are expected to be larger. Fortunately, we were 
able to detect both 1665 and 1667 MHz OH masers. As was 
found in vLOO, the most blue-shifted maser features in both tran- 
sititions are coincident, and correspond to the amplified stel- 
lar image. We thus estimate an additional factor in the posi- 
tional uncertainty from the internal scatter between the 1665 
and 1667 MHz most blue-shifted maser features. This additional 
scatter factor, of the order of ~ 0.6 mas, was also added quadrati- 
cally to the positional error estimate described above, producing 
a typical error of ~ 0.7 mas. As a result, the;t/^-fit to both 1667 
and 1665 MHz OH masers of U Her simultaneously, has a re- 
duced - 1.61 for 13 degrees of freedom. These results 



imply that the error estimates for all three of the sources are a 
fairly good representation of the true astrometric errors. 

Finally, we have checked the robustness of the astrometric 
results by randomly removing individual data points and per- 
forming new fits on the reduced data set. We find that for S CrB 
and RR Aql the fit solutions are consistent within Icr. For U Her 
the solutions are consistent within ~ 1.2cr, hinting at a slight 
underestimate of the fitting errors for U Her. 

3. Results 

The results of the;(f^-fits are presented in Table|2] The table con- 
tains the parallax and proper motion results derived from the re- 
cent observational epochs as well as, for comparison, the results 
from V03 for U Her and S CrB. For these sources we have also 
performed a fit combining the entire data-set, allowing for an ad- 
ditional zero-position off-set between the earlier epochs and the 
data presented here. Such an additional off-set is needed because 
of changes in the calibrator situation and since, especially for 
S CrB, we might be tracing a new maser feature. Additionally, 
the errors for the early epochs in the combined fitting have been 
increased by a factor of two, as in V03 only the formal position 
fitting errors were used while in the original scheme the iono- 
spheric contributions could not be easily estimated easily esti- 
mated. The table gives the velocity Vlsr of the maser feature 
used for the fit with respect to the Local Standard of Rest (LSR) 
as well as which maser transition (1665 or 1667 MHz) was used 
for the fit. 

3.1. RRAqI 

RR Aql is one of the new sources added to the sample, for which 
Hipparcos did not manage to obtain a significant parallax. We 
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Table 2. Results 



Source 


Maser 
(MHz) 


Note" 


VLBI parallax 
(mas) 


VLBI Proper motion 
RA, dec (mas/yr) 


# epochs 


I'LSR 

(km s-') 


<^VLBI 

(pc) 


SCrB 


1667 
both 
both 


A 
B 


2.39 + 0.17 
2.31+0.33 
2.36 + 0.23 


-8.58 + 0.38 -13.21 +0.61 
-9.08 + 0.27 -12.49 + 0.33 
-9.06 + 0.23 -12.52 + 0.29 


5 
8 
5+8 


4.6 km s"' 
3.2/2.9 km S-' 


418!- 


UHer 


both 
1667 
both 


A 
B 


3.76 + 0.27 
3.61 + 1.04 
3.74 + 0.61 


-16.99 + 0.77 -11.88 + 0.50 
-14.94 + 0.38 -9.17 + 0.42 
-14.98 + 0.29 -9.23 + 0.32 


5 
12 
5+12 


-20.2 km S-' 
-20.4 km S-' 


266!^^ 


RR Aql 


1665 




1.58 + 0.40 


-25.11 +0.74 -49.82 + 0.54 


5 


24 km s"' 


633_|28 



" A: results from V03; B: results from combined fit including V03 epochs 
* LSR Velocity of fitted maser features 
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Fig. 4. Period vs. Ko magnitudes for the OH maser stars with 
VLBI interferometric distances, including those from V03. The 
soUd symbols and the error bars are determined using the VLBI 
distances. The open symbols are the values using the Hipparcos 
distances. The solid line is the P-L relation determined by 
IWhitelock & Feasli ( 120001) on the oxygen rich Mira stars ob- 
served with Hipparcos, the dashed lines are the spread in the 
relation due to the error in the P-L relation zero-point. 



did not manage to detect any 1667 MHz maser emission, but 
were fortunately able to image a ~ 0. 1 Jy/beam 1665 MHz maser 
feature in the low spectral resolution data. The feature was found 
at 24 km s"\ blue-shifted with respect to the stellar velocity, but 
at a velocity unlikely to belong to an amplified stellar image. We 
managed to detect and image the feature at 5 epochs with respect 
to the nearby in-beam calibrator. This yielded the smallest par- 
allax determined to date for OH maser astrometry, even though 
the errors are dominated by the fact that the maser is weak and 
is only detected in the low resolution spectral line data-set. 

3.2. S CrB 

S CrB was previously monitored with the VLBA (V03), at 
which time we were able to detect bright masers at both 1665 



and 1667 MHz (at Vlsr = 3.2 and 2.9 km s^' respectively). 
In the new epochs we were only able to detect a 1667 MHz 
maser feature, red-shifted with respect to the stellar velocity, at 
Vlsr - 4.6 km s"' . The feature was fairly weak during the first 2 
epochs, with the flux decreasing from 0.07 Jy/beam in the first to 
below 0.05 Jy/beam in the second epoch. The third epoch turned 
out to be unusable, but during the last three epochs the maser 
had brightened again to ~ 0.11 Jy/beam. Thus, the astromet- 
ric errors were almost a factor of two smaller in the last three 
epochs. Both the parallax and proper motion fit results from the 
newest epochs are fully consistent with the result obtained from 
the earlier epochs. 

3.3. U Her 

U Her is the OH maser source which has been the target of the 
longest astrometric VLBI monitoring campaign, and cuiTently 
has data spanning almost 12 years (vLOO, V03). In the new ob- 
servations, we again detected the most blue-shifted maser fea- 
ture corresponding to the amplified stellar image in both the 
1665 and 1667 MHz maser transition, with a flux varying be- 
tween 0.8 and 1 . 1 Jy/beam. In addition, we detected a number of 
weaker features which are discussed in § 14.21 While the paral- 
lax is also fully consistent with the earlier determined parallax, 
albeit with a significantly smaller error, the proper motion has 
changed considerably. As was already found in V03, the cah- 
brator used for phase referencing, J1630-I-2131, shows structural 
changes on timescales at the order of several months. This is 
possibly the cause of the differences in proper motions. FigfTOlof 
the online material shows a series of images of Jl 630+2131 for 
each observational epoch indicating these changes. As described 
in V03, the internal calibrator motions are expected to be rela- 
tively slow compared to the annual timescale of the parallax and 
we thus think it unlikely that this will bias the parallax measure- 
ment. However, the motions will introduce additional scatter and 
some caution has to be taken with respect to the quoted uncer- 
tainties. Because of the proper motion difference, caution has to 
be taken in using the fitted motion resulting from the combined 
epochs fitting. 

4. Discussion 

4.1. Comparisson with previous results 

The distances to the OH bearing Mira variable stars S CrB and 
U Her derived in this paper are fully consistent with the previ- 
ously published values in V03. The uncertainty on the parallax 
measurement decreased by over a factor of 2, even with signif- 
icantly fewer observational epochs, illustrating the astrometric 
improvements. The eiTors on the proper motion measurements 
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are larger, mainly due to the shorter time-baseline of the new 
observations. The most accurate proper motions have therefore 
been obtained by simultaneously fitting the maser positions from 
both the current and previous observational epochs as described 
above. 

The astrometric measurements are also a significant im- 
provement on the Hipparcos derived values, which suffered from 
the fact that the stars in our sample are relatively faint, vari- 
able, AGB stars (see V03 for more details). When comparing the 
VLBI derived distances in Table|2]with distances derived from a 
P-L relation in Table [1] we find that, while for both S CrB and 
RR Aql the distances are reasonably consistent, the U Her VLBI 
distance is significantly lower. This can also be seen in Fig. |4] 
which shows position of the OH maser stars with VLBI paral- 
laxes in the P-L diagram. As was speculated in vLOO and V03, 
the discrepancy for U Her could be a result from U Her pulsating 
in the fundamental mode, while the other stars are pulsating in 
the first overtone mode. However, as it has been argued that al l 
Mira stars are pulsating in the fundamental mode (IWoodl[T998l) . 
the offset of U Her more likely highlights the scatter in the P-L 
relation or the uncertainty in individual magnitude determina- 
tions. We rule out that the offset is introduced by the absence of 
an in-beam calibrator for this source, because we can imagine nu 
systematic effect that can bias the parallax to larger values over 
all the years we have monitored the source. 

4.2. Morphology and internal motions of the OH masers of 
UHer 

While for S CrB and RR Aql we only detect maser features in 
a few neighbouring red-shifted and blue-shifted spectral chan- 
nels respectively, we observe a large number of both blue- and 
red-shifted features for U Her. Fig. |3a-e) show the first mo- 
ment greyscale maps with overplotted total intensity contours 
for the 5 succesful U Her epochs. The maps are all referenced 
the most blue-shifted maser feature that has been identified as 
the stellar image. Several of the maser features are detected in 
multiple epochs and these are shown in Fig. 13 f). The brightest 
red-shifted masers are all detected with Vlsr between -10 and 
-12 km s~^ while the blue-shifted masers are located at Vlsr be- 
tween - 1 7 and -2 1 km s" ' . There is no obvious difference in the 
extent of the red- and blue-shifted maser features, which are both 
spread over ~ 200 x 200 mas, corresponding to ~ 55 x 55 AU at 
the 266 pc distance of U Her Fig.|5]does show that the maser dis- 
tribution is distinctly asymmetrical with respect to the stellar im- 
age, with hardly any maser features detected to the North-East. 
This is significantly different from the 1665 and 1 667 GHz maser 
distribution observed in 1984 with MERLIN (Chapma n et al.L 
11994) . but similar to our previous 1665 and 1667 GHZ OH 
maser observations between 1995 and 2002 (vLOO, V03). As the 
22 GHz H2O maser emission is also clearly weakest towards the 
North-East ( Vlemming s et all [2002: Bains et al., 2003), the ob- 
served non-spherical maser distribution is likely due to a real 
asymmetry in the CSE. 

As several U Her maser features, besides the stellar image, 
can be traced over multiple epochs, we are able to study the 
internal motions in the maser envelope. One red-shifted maser 
feature, at Vlsr - -11.51 km s"', was persistent over the en- 
tire year of observation, while three maser features are found in 
4 epochs. A further 5 features can be traced over three epochs 
while another 5 are detected in only 2 epochs. Although none 
of these features are bright enough to be included in the regu- 
lar astrometric fitting, we are able to test the consistency of our 
observations by determining the motions of these features after 



self-calibrating on the blue-shifted maser feature amplifying the 
stellar emission. By applying the self-calibration we remove the 
stellar proper motion and parallax, after which we fit for resid- 
ual proper motions of the maser features. For most of the maser 
features we do not detect significant proper motion. Only for 
one red-shifted maser feature, that was detected in the first three 
epochs at Vlsr = -10.81 km s"', we detect a proper motion 
of fi - 5.1 + 1.5 mas yr"\ which corresponds to a velocity of 
6.4 +1.9 km s^'at the distance to U Her. Assuming a spheri- 
cally expanding shell, this implies, from the angular off-set to 
the stellar image and the radial velocity, an outflow velocity of 
~ 7.5 km s 'in the U Her OH maser region. The proper motion 
of this feature is indicated with a thick line in Fig.|5] Performing 
a weighted ensemble average of the masers in three quadrants 
(ff < and 6 > 0, a < Q and 6 < 0, a > and 6 < 0) reveals no 
significant motions. Considering that we only observe the front 
and back side of the maser shell this is not surprising, as the out- 
flow will be predominantly in the radial direction. Thus, as no 
residual proper motion is detected, the motion determined from 
the astrometric fitting to the amplified stellar image can be taken 
as the true stellar proper motion, with only a possible contribu- 
tion from the internal motions of the phase reference source. 

As we have concluded that the majority of the maser spots 
show no significant proper motion when referencing them to 
the stellar image, the residual motion of the maser features is 
likely due to a combination of the intrinsic astrometric errors 
and turbulent motions in the masing medium. We find that the 
rms scatter of the masers detected at multiple epochs is o-pos = 
2.0 ±1.3 mas. Taking into account the intrinsic positional fit- 
ting uncertainty, including the error discussed in ^2.41 which for 
these weak maser features is of the order of 1 .5 mas, small but 
significant residual motions are present . These we attribute to 
turbulence at the level of 1 - 2 km s only sli ghtly higher t han 
the typical assumed value of 1 km s~' (Diamond et aUll985h . 

5. Concluding remarks 

We have been able to significantly improve the parallax mea- 
surements for a number of OH maser stars. The improvements 
were the results of using nearby phase reference sources in the 
form of in-beam calibrators, and generally good observing con- 
ditions due to a minimum in solar activity. As a result the paral- 
lax uncertainties were decreased by up to a factor of ~ 4 to the 
order of ~ 0.2 mas, using only 5 observing epochs compared to 
the 8-12 epochs presented in V03. Thus, we have convincingly 
shown the capability of VLBI astrometry to determine distances 
up to 1 kpc to AGB stars using their circumstellar OH masers, 
even when no special maser feature amplifying the stellar image 
is found. The limit on distance measurements could potentially 
be stretched to ~ 2 kpc or :7r ~ 0.5 mas when the OH maser is 
strong (> 0.25 Jy/beam) and a bright compact in-beam calibra- 
tor is available. The precision will be further increased when a 
number of different maser features can be traced over a sufficient 
number of epochs. This opens up the possibility to use OH maser 
astrometry to determine the parallaxes to a large number of OH 
maser stars with previously unknown distances, which can for 
instance be used to contruct a P-L relation for the most heavily 
enshrouded AGB stars. 
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Fig. 5. (a-e) Images for 5 epochs of the 1667 MHz OH maser emission of U Her referenced to the most brightest maser. The 
brightest and most blue-shifted maser feature is prominent in all images and is identified with the stellar image amplified by the 
masing material at the front side of the expanding shell (V03). The greyscale indicates the first moment map with the velocity scale 
indicated at the top of the figures. The contours indicate the total intensity I, drawn at levels of 40, 80, 160 and 320 mJy/Beam. The 
rms noise varies somewhat for the difl'erent epochs but is ~ 12 mJy/Beam. The peak intensity is given for each individual image, 
(f) The maser features that were observed during at least 2 epochs with vectors indicating their proper motion. The arrow denotes a 
3cr proper motion detection, no significant proper motion was detected for any of the other features but for illustration purposes we 
show l-3cr proper motions with the regular vectors. The solid square indicates the position of the stellar image. 



the mixed bandwidth calibration scheme in ParselTongue. 
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Fig. 6. Results from the 8.4 GHz VLA in-beam calibrator search for RR Aql and S CrB. The images are labelled according to 
Tableland include non-detections for completeness. No special calibration to improve the images was performed. Peak flux and 
contour levels are given for each image individually. For RR Aql and S CrB, VLBA in-beam calibration was sucessfully performed 
on the sources RRAQL.03 and SCRB.02. 
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Fig. 7. As Fig. |6] for U Her and R Aql. VLBA in-beam calibration was unsuccesfuUy attempted on UHER.02. No masers were 
detected for R Aql during the VLBA observations. 
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Fig. 8. As Fig.|6]for R Cas and Y Cas. No masers were detected during the VLBA observations for eitlier star. 
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Fig. 9. As Fig.[6]for RT Oph and W Hya. None of the sources were of sufficient quality to be used as in-beam calibrators. 
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Table 3. X-Band VLA observation results: R esults from the X-ba nd VLA search for suitable in-beam calibrators. The sources 
were picked from the 1.4 GHz NVSS survey dCondon et all [1998 ^. The observations were taken in D-configuration at May 21st 
2003. The table gives the name identifying the nearby maser source with the NVSS name and integrated flux. Additionally, we give 
the coordinates of the peak emission, and the integrated and peak 8.4 GHz flux determined using the AIPS task JMFIT, when a 
single component Gaussian fit was possible for the strongest component. The last column contains the seperation Afl between the 
continuum source and the maser star in arcminutes. For those sources that were used as in-beam calibrators and detected in the 
VLBA L-band observations we also give 7MF/r intergrated fluxes. 



Source 


NVSS name 


NVSS Flux 
(mJy) 


hh mm ss 


Snooo 


Flux 
(mJy) 


Peak Flux 
(mjy/beam) 


^e 

(arcmin) 


RRAQL.Ol 
RRAQL.02 
RRAQL.OS" 

RRAQL.04 
RRAQL.05 


NVSS J195826-020352 
NVSS J195759-014012 
NVSS J195655-013615 
L-band VLBI 
NVSS J195728-021337 
NVSS J195627-014425 


116.8 ±3.5 
115.5 ±4.1 
94.3 ± 2.9 
48.1 ±0.2 
36.6 ± 1.2 
16.0 ±0.6 


19 58 26.5826 
19 57 59.5262 
19 56 55.2046 

19 56 27.0759 


-02 03 51.833 
-01 40 14.670 
-01 36 15.426 

-01 44 27.040 


13.57 ±0.19 
21.00 ±0.08 
70.96 ±0.12 

7.20 ± 0.04 


12.77 ±0.11 
20.48 ± 0.05 
69.95 ± 0.07 

7.14 ±0.02 


16.5 
14.2 
19.8 

20.5 
19.3 


SCRB.Ol 
SCRB.02°'' 

SCRB.03 
SCRB.04 
SCRB.05 
SCRB.06 


NVSS J152101+311537 
NVSS J152209+314414 
L-band VLBI 
NVSS J152233+310911 
NVSS J152109+310334 
NVSS J152012+310738 
NVSS J152111+312816 


175.5 ± 5.3 
347.2 ± 10.4 
252.7 ±0.1 
51.3 ± 1.9 
21.3 ±0.7 
11.9 ± 0.5 
11.8±0.5 


15 21 01.2872 
15 22 09.9937 

15 21 09.7814 
15 20 12.7797 


+31 15 37.796 
+31 44 14.361 

+31 03 34.369 
+31 07 37.135 


97.48 ± 0.69 
465.18 ±2.95 

2.44 ±0.11 
1.27 ±0.03 


96.99 ± 0.40 
458.56 ± 1.69 

1.62 ±0.05 
1.11 ±0.01 


8.1 
24.2 

19.7 
18.7 
20.9 
6.8 


UHER.Ol 
UHER.02'' 

UHER.03 
UHER.04 
UHER.05 
UHER.06 
UHER.07 


NVSS J162605+191624 
NVSS J162539+183727 
L-band VLBI 
NVSS J 162435 +185 344 
NVSS J 162544+ 184406 
NVSS J162508+191647 
NVSS J162508+190212 
NVSS J162650+191246 


1 14.0 ± 3.4 

32.8 ± 1.1 
4.5 ±0.3 

29.9 ± 1.0 
20.2 ± 0.7 
20.0 ± 1.1 
7.8 ±0.5 
12.8 ±0.6 


16 26 05.5769 
16 25 38.9655 

16 24 35.6826 
16 26 50.0838 


+ 19 16 24.999 
+ 18 37 25.691 

+ 18 53 45.245 
+ 16 26 50.0838 


9.24 ± 0.07 
5.10 ±0.02 

0.78 ± 0.10 
0.94 ± 0.06 


2.63 ± 0.02 
3.92 ±0.01 

0.29 ± 0.03 
0.37 ± 0.02 


23.2 
16.2 

17.0 
9.5 
25.0 
12.7 
24.3 


RAQL.Ol" 

RAQL.02 
RAQL.03 
RAQL.04 


NVSS J 1907 15 +075240 
L-band VLBI 
NVSS J190507+082958 
NVSS J190552+075754 
NVSS J190502+080239 


268.1 ± 9.5 
1.1 ±0.2 
200.0 ± 6.0 
105.6 ± 3.2 
56.0 ± 3.5 


19 07 15.5506 

19 05 07.0343 
19 05 52.4491 


+07 52 43.409 

+08 29 57.507 
+07 57 53.358 


35.07 ± 0.28 

8.95 ± 0.07 
6.88 ± 0.05 


20.10 ±0.11 

7.66 ± 0.04 
3.80 ± 0.02 


24.9 

24.6 
17.5 
22.7 


RCAS.Ol 

RCAS.02 

RCAS.03 

RCAS.04'' 

RCAS.05 


NVSS J235918+513815 
NVSS J235959+5 11727 
NVSS J235908+5 10703 
NVSS J000019+5 11254 
NVSS J235735+5 10122 


58.9 ± 1.8 

8.7 ±0.5 
7.3 ±0.5 
6.9 ±0.5 

6.8 ±0.5 


23 59 19.0708 

23 59 09.8053 
00 00 19.3222 
23 57 34.8546 


+51 38 15.198 

+51 07 08.929 
+51 12 54.106 
+51 01 20.852 


9.97 ± 0.07 

0.50 ± 0.04 
17.49 ±0.17 
10.03 ± 0.06 


8.11 ±0.03 

0.47 ± 0.02 
17.60 ± 0.10 
10.06 ± 0.03 


17.1 
15.9 
17.7 
20.7 
23.3 


YCAS.Ol" 

YCAS.02 
YCAS.03 
YCAS.04 


NVSS J000108+552638 
L-band VLBI 
NVSS JOOO 156+554730 
NVSS J000344+555625 
NVSS J000339+552421 


289.7 ± 8.7 

68.5 ± 0.3 

115.8 ±4.2 
84.4 ± 2.6 

42.6 ± 1.7 


00 01 08.7229 

00 01 56.5666 
00 03 44.1589 


+55 26 38.668 

+55 47 26.041 
+55 56 26.483 


33.93 ±0.11 

24.82 ±0.10 
7.80 ± 0.07 


31.76 ± 0.06 

24.56 ± 0.06 
2.68 ± 0.02 


23.5 

13.7 
15.9 
16.7 


RTOPH.Ol 
RTOPH.02 
RTOPH.03 
RTOPH.04 


NVSS J175757+111318 
NVSS J175637+104729 
NVSS J175630+1 11433 
NVSS J175613+1 10835 


58.8 ±2.1 

37.7 ± 1.2 

22.8 ±0.8 
16.1 ±0.6 


17 56 37.9921 
17 56 30.5689 


+ 10 47 28.239 
+ 11 14 33.598 


3.15 ±0.05 
3.77 ± 0.03 


2.70 ± 0.03 
3.70 ± 0.02 


21.1 
22.7 
4.4 
4.9 


WHYA.Ol 
WHYA.02 
WHYA.03 
WHYA.04 
WHYA.05 


NVSS J135030-281027 
NVSS J134845-281239 
NVSS J134740-282942 
NVSS J134832-283356 
NVSS J134833-280909 


16.3 ±0.7 
9.9 ±0.6 
8.4 ± 0.5 
7.7 ±0.5 
7.3 ±0.5 


13 50 30.2936 
13 48 44.9729 

13 48 32.2336 
13 48 33.6641 


-28 10 27.173 
-28 12 42.180 

-28 33 56.579 
-28 09 08.811 


2.97 ± 0.09 
1.49 ±0.05 

1.08 ±0.03 
0.67 ± 0.04 


2.56 ± 0.05 
1.41 ±0.03 

1.08 ±0.02 
0.54 ± 0.02 


22.6 
10.1 
19.6 
13.6 
14.3 



° In-beam calibrator used during VLBA observations. 

* Also known as JVAS source J1522+3144 with previously know VLBI position. 
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Fig. 10. The primary phase reference calibrator J 1 630+21 3 1 used for U Her. The panels are labelled by the number of days after the 
first epoch of observations and include epoch e, for which we were unable to properly image the OH masers. The contours start at 
1 mJy and each subsequent contour level is multiplied by a factor of "\^. 



